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ABSTRACT 

Numerical simulations have shown that the often cited radiation pressure barrier to accretion 
onto massive stars can be circumvented, when the radiation field is highly anisotropic in the 
presence of a circumstellar accretion disk with high optical depth. Here, these studies of the 
so-called flashlight effect are expanded by including the opacity of the innermost dust-free but 
potentially optically thick gas regions around forming massive stars. In addition to frequency- 
dependent opacities for the dust grains, we use temperature- and density-dependent Planck- and 
Rosseland mean opacities for the gas. The simulations show that the innermost dust-free parts 
of the accretion disks are optically thick to the stellar radiation over a substantial fraction of the 
solid angle above and below the disk's midplane. The temperature in the shielded disk region 
decreases faster with radius than in a comparison simulation with a lower constant gas opacity, 
and the dust sublimation front is shifted to smaller radii. The shielding by the dust-free gas 
in the inner disk thus contributes to an enhanced flashlight effect, which ultimately results in a 
smaller opening angle of the radiation pressure driven outflow and in a much longer timescale of 
sustained feeding of the circumstellar disk by the molecular cloud core. We conclude that it is 
necessary to properly account for the opacity of the inner dust-free disk regions around forming 
massive stars in order to correctly assess the effectiveness of the flashlight effect, the opening 
angle of radiation pressure driven outflows, and the lifetime and morphological evolution of the 
accretion disk. 

Subject headings: ISM: dust, extinction — ISM: molecules — methods: numerical — stars: circumstellar matter — 
stars: formation — stars: winds, outflows 
©2012. All rights reserved 



1. Introduction 

During their lifetime massive stars can exert a radia- 
tive force on their surroundings, which can be higher 
than the gravitational attraction. How these massive 
stars can sustain accretion of high-opacity dusty gas 
was an open question for decades. Semi-analytical 
(Kahn 1974) and first radiation-hydrodynamical com- 
putations in spherical symmetry (Yorke & Krugel 
1977) supported the idea that the radiative force of 



a star slows and could ultimately prevent the accre- 
tion of material from the molecular core, leading to 
an upper stellar mass limit of < 40 M Q . Wolfire & 
Cassinelli (1987) showed how the interstellar medium 
dust mixture had to be altered to allow the formation of 
the most massive stars known by spherical symmetric 
accretion flows. They found that the absorption prob- 
ability of the stellar environment has to be reduced to 
unlikely low values, e.g. by reducing the dust to gas 
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mass ratio significantly and removing large graphite 
grains. 

Based on the fact that in all of these one-dimensional 
simulations the collapse was reversed by the re-emitted 
thermal dust radiation, Nakano (1989) inferred that the 
radiative force onto the accretion flow can be dimin- 
ished by an anisotropic optical depth of the environ- 
ment. Such an anisotropy can naturally be produced 
by a circumstellar accretion disk, forming during the 
pre-stellar core collapse due to angular momentum 
conservation. 

The first numerical study aimed at verifying this 
idea was by Yorke & Sonnhalter (2002), who further 
expanded on the concept of the so-called "flashlight 
effect" for the anisotropic thermal dust radiation, as 
introduced by Yorke & Bodenheimer (1999). Whereas 
Yorke & Bodenheimer (1999) could verify the signifi- 
cance of the flashlight effect, in particular for estimates 
of luminosity and spectral appearance, they only con- 
sidered molecular core masses less than 10 M G . In the 
Yorke & Sonnhalter (2002) simulations higher initial 
core masses up to 120 M were considered, but the op- 
tical depth of the forming accretion disk turned out to 
be too low to provide a sufficient anisotropy of the radi- 
ation field. As a result, the accretion onto the protostar 
stopped shortly after the disk formation, limiting the 
maximum final stellar mass to 42.9 M . In Krumholz 
et al. (2009), the authors claimed further feeding of the 
star and its circumstellar disk by a radiative Rayleigh- 
Taylor instability in the bipolar cavity regions. The 
most massive (and still accreting) protostar formed in 
this simulation was 41.5 M G . 

All in all, no star much more massive than the 
ID radiation pressure barrier of M™ ax m 40 M was 
formed in these numerical simulations. 

In Kuiper et al. (2012) the radiation transport 
method used in the Krumholz et al. (2009) study 
(gray flux-limited diffusion = FLD approximation) 
was compared with a much more sophisticated hybrid 
radiation transport approach, namely a combination of 
a frequency-dependent ray-tracing scheme for the stel- 
lar irradiation and a gray FLD solver for the thermal 
dust (re-)emission. Kuiper et al. (2012) showed that 
the radiation-pressure-dominated cavities became ra- 
diative Rayleigh-Taylor unstable only in the runs using 
the gray FLD approximation for the stellar feedback, 
whereas they remained stable in the runs using the 
hybrid approach. These results were backed up by an- 
alytical estimates of the radiative forces acting on the 
expanding cavity shell. These results suggest that a ra- 



diative Rayleigh-Taylor instability is not a mechanism 
for overcoming the radiation pressure barrier. 

The next numerical study aimed at verifying the 
idea of Nakano (1989) regarding the anisotropic radia- 
tion field effect was presented by Kuiper et al. (2010a). 
The radiation hydrodynamical simulations presented 
therein confirmed the radiation pressure barrier for 
spherically symmetric flows and demonstrated that this 
problem can be circumvented by an anisotropic ther- 
mal radiation field, also known as the flashlight effect, 
as proposed. It was shown that the flashlight effect 
even allows the most massive stars known to be formed 
by disk accretion despite of their strong radiation feed- 
back. Furthermore, because sufficient angular momen- 
tum transport can be provided by gravitational torques 
forming in the self-gravitating accretion disk (Kuiper 
et al. 201 1), there is no angular momentum barrier to 
accretion onto the protostar. 

Thus, circumstellar accretion disks appear to be a 
necessary ingredient for massive star formation. For 
observational techniques and first disk candidates, we 
refer the reader to the reviews by Zhang (2005) on ob- 
servations and Cesaroni et al. (2007) on theory and ob- 
servations of disks around massive stars. 

Due to the fact that the optical depth of the dusty 
regions around massive stars will be highly dominated 
by dust opacities rather than gas opacities, the numer- 
ical studies mentioned above have used simplified gas 
opacities. In Yorke & Sonnhalter (2002), the dust-free 
gaseous regions around the forming massive star were 
- especially in the disk's midplane - most likely hid- 
den in the central sink cell of the size of 40, 80, and 
160 AU for an initial core mass of 30, 60, and 120 M . 
The sink cells were treated as being optically thin. 
Krumholz et al. (2009) used a constant Planck- and 
Rosseland mean gas opacity of Kp — kr = 0.1 cm 2 g" 1 
and a resolution of the highest refinement level of the 
AMR grid of 10 AU, i.e. not resolving the scale height 
of the dust-free gas disks. Kuiper et al. (2010a, 2011, 
2012) used a constant Planck- and Rosseland mean gas 
opacity of Kj> = acr = 0.01 cm 2 g -1 and the resolution 
of the inner gas disk was about 1 AU at the inner com- 
putational boundary, decreasing to larger radii propor- 
tional to the radius. 

Here, we have improved on this detail by newly 
implementing the temperature- and density-dependent 
Planck- and Rosseland mean gas opacities of Helling 
et al. (2000). For the dust grains we use the frequency- 
dependent opacities by Laor & Draine (1993), which 
were used in our previous simulations. In addition to 
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our default numerical grid of the previous simulations, 
we here also show the results of extending the compu- 
tational domain towards a smaller size of the central 
sink cell of R m [ n = 5 AU, including parts of the dust- 
free accretion disk around the forming massive star 
and increasing the polar resolution at the inner bound- 
ary to 0.5 AU. 

To our best knowledge, there are no previous dy- 
namical studies of the gaseous disks around massive 
stars in the literature. In Tanaka & Nakamoto (201 1), 
the authors investigated in analytical and numerical 
calculations the radiative force by the direct stellar ir- 
radiation of the dust sublimation front and in the in- 
ner dust-free gas disk for a given static disk configu- 
ration. They propose that the dust sublimation front in 
the disk's midplane will be shielded by the optically 
thick inner gas disk. 

Here, we study this scenario in dynamical simula- 
tions of the formation and evolution of circumstellar 
accretion disks. In addition to the direct stellar irradi- 
ation we include the thermal (re-)emission of both the 
dusty gas and the dust-free gas and investigate the dy- 
namical consequences of the optical depth of the dust- 
free regions for the dust sublimation front, the accre- 
tion flow, the flashlight effect, and the bipolar outflow 
regions. In high-mass star forming regions, the radia- 
tive force potentially causes or influences the observed 
wide angle outflows, although it remains unclear, if 
these outflows are mainly driven by radiation or by a 
combination of magnetic fields and centrifugal forces. 

In Sect. 2 we describe the changes made to the nu- 
merical code with an emphasis on the newly imple- 
mented gas opacities by Helling et al. (2000). Sect. 3 
provides an overview of the cases considered, their ini- 
tial conditions, and the numerical configuration. We 
present the results of these simulations in Sect. 4, giv- 
ing the physical explanations and discussion of the 
consequences in Sect. 5. We close this report in Sect. 7 
with a discussion of the limitations to the current study 
and an outlook for future investigations. 

2. Methods 

The numerical code used to run the pre-stellar 
core collapse simulations is based on the one used 
in previous studies (Kuiper et al. 2010a, 2011, 2012). 
For hydrodynamics we use the open source magneto- 
hydrodynamics code Pluto (Mignone et al. 2007). For 
a detailed description of the derivation, numerical im- 
plementation, and benchmarking of the radiation trans- 
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Fig. 1. — Planck mean gas opacities by Helling et al. 
(2000) as a function of the gas temperature for two 
values of the gas density, p = 10~ 11 gem -3 (solid 
line) and p = 10~ 20 g cm 3 (dashed line). The hor- 
izontal line marks the constant gas opacity of k — 
0.01 cm 2 g -1 used in the comparison simulation run. 
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Fig. 2. — Rosseland mean gas opacities by Helling 
et al. (2000) as a function of the gas temperature for 
two values of the gas density, p = 10~ n g cirT 3 (solid 
line) and p = 10~ 20 g cirT 3 (dashed line). The hor- 
izontal line marks the constant gas opacity of k — 
0.01 cm 2 g -1 used in the comparison simulation run. 
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port solver in use, please see Kuiper et al. (2010b). The 
evolution of the central star is computed via fits to the 
stellar evolutionary tracks of Hosokawa & Omukai 
(2009) including the effect of accretion onto forming 
high-mass stars. As in our previous studies, the opac- 
ity of dust grains is computed taking into account the 
frequency-dependence for the stellar irradiation and 
using frequency-averaged (gray) Planck- and Rosse- 
land mean opacities for the thermal dust emission. We 
use the frequency-dependent dust opacities by Laor & 
Draine (1993). 

Here, we improved the radiation transport by im- 
plementing gray, temperature- and density-dependent 
Planck- and Rosseland mean opacities of the gas. 
Whereas dust grains will dominate the opacity in ma- 
jor parts of the collapsing pre-stellar core, gas opacities 
could become important in the dust-free regions close 
to the protostar. Here, we have implemented the gas 
opacities computed in Helling et al. (2000). These gas 
opacities were also used in a study deriving opacities 
for a broad parameter range (in density and temper- 
ature) for proto-planetary disks in general (Semenov 
et al. 2003). Opacity values for given logarithmic gas 
temperatures and densities are interpolated bi-linearly 
from a 71 x 71 grid. The Planck- and Rosseland mean 
gas opacities as a function of gas temperature are de- 
picted in Figs. 1 and 2 for two values of gas density. 
The optical depth of regions of lower temperatures 
than plotted will most likely be dominated by dust 
grains. 

3. Simulations overview, initial conditions, and 
numerical configuration 

We have performed three different simulations: 
For the "fiducial run" we computed the collapse of a 
pre-stellar core including the newly implemented gas 
opacities. A "comparison run" with the same initial 
conditions and numerical configuration is performed 
using a constant gas opacity of k — 0.01 cm 2 g _1 . 
To determine the dependence on the extent of the 
computational domain towards the direction of the 
central star, the fiducial run is repeated (called the 
"convergence run") using a smaller sink cell size of 
^min = 5 AU instead of the value R m \ n = 10 AU used 
in the fiducial and comparison runs. An overview of 
these simulations is given in Table 1 . 

The initial condition of the pre-stellar cores is given 
by an outer core radius of 0.1 pc, a mass of 100 M , 
and a temperature of 20 K. The core is initially in rigid 



solid body rotation. 

The numerical grid in use is identical to several of 
our previous studies. The resolution of the spherical 
grid at the inner computational boundary is about 1 AU 
in the radial as well as in the polar direction, assuming 
axial symmetry. The polar extent of the grid covers the 
angle from 0° to 90°, assuming midplane symmetry. 
Towards the outer computational boundary, the grid 
cell size increases linearly with the radius. In the run 
using the smaller central sink cell size of R m \ n = 5 AU, 
the radial grid is extended towards the smaller radii 
using a uniform step size of 1 AU. For more details 
on the numerical grid, the numerical solvers, and the 
sub-grid models such as shear viscosity and the dust 
model, please see Kuiper et al. (2010b) and Kuiper 
etal. (2010a). 



Run tag 


tfmin [AU] 


""gas 


M. [M Q ] 


fiducial run 


10 


H 


60.5+ 


comparison run 


10 


c 


52.1 


convergence run 


5 


H 


58.0+ 



Table 1: Overview of simulations performed. From left 
to right the columns denote the name tag, the radius of 
the inner boundary condition, and the gas opacities in 
use: "H" refers to the usage of the frequency-averaged 
Planck- and Rosseland mean gas opacities by Helling 
et al. (2000), "c" refers to the usage of a constant gas 
opacity of 0.01 cm 2 g . The last column denotes the 
resulting final mass of the star forming in the central 
sink cell. For the fiducial and convergence run the 
"+" sign denotes the fact that in the end of the simula- 
tion the star is still accreting. The remaining mass in 
the computational domain in the fiducial run is 4.0 M 
at t = 253 kyr. The remaining mass in the compu- 
tational domain in the convergence run is 20.8 M at 
t = 80 kyr. 

4. Results 

During the initial free-fall epoch and the early 
epoch of disk formation and evolution, the fiducial, 
comparison, and convergence runs show no qualitative 
and only minor quantitative differences. The central 
star eventually becomes so luminous that its radiative 
force exerted onto the environment overcomes grav- 
ity. Over time, these radiative forces lead to signif- 
icant mass loss from the pre-stellar molecular core. 
The evolution of the simulation runs using differ- 
ent gas opacities starts to deviate significantly once 
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the radiation-pressure-dominated cavities are formed. 
First, the opening angle of the radiation-driven outflow 
is smaller for the run with the newly implemented gas 
opacities, resulting in different mass loss rates from 
the pre-stellar core, see Fig. 3. This mass loss rate 
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Fig. 3. — Mass loss of the pre-stellar molecular cores 
due to radiative acceleration as a function of time. The 
solid line denotes the fiducial simulation run with gas 
opacities by Helling et al. (2000). The dashed lines de- 
notes the comparison run using a constant gas opacity. 

is computed by integrating the mass flux through the 
outer computational boundary at R max = 0.1 pc. This 
boundary is treated as a semi-permeable wall, i.e. ma- 
terial with outward radial velocity is allowed to leave 
the computational domain, but inflow of material is 
prohibited. The fiducial run with the newly imple- 
mented gas opacities initially has a much lower mass 
loss rate. Second, the fiducial and the comparison run 
have very different accretion histories of the forming 
star, see Fig. 4. The fiducial run has a much longer 
epoch of sustained accretion than in the comparison 
run (see explanation in Sect. 5). 

Radiative acceleration stops disk accretion in the 
comparison run at about 65 kyr. After this initial shut 
off the mass in the remnant disk piles up and leads to 
an epoch of low accretion (a few 10~ 5 M Q yr _1 ) from 
roughly 70 to 90 kyr. Eventually, the radiation force is 
able to disperse the remnant disk and accretion shuts 
off. 
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Fig. 4. — Mass of the central star (upper panel) and 
accretion rate onto the star (lower panel) as a function 
of time. The solid lines denote the fiducial simulation 
run with the gas opacities by Helling et al. (2000). The 
dashed lines denote the comparison run using a con- 
stant gas opacity. 
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5. Discussion 

In the following, we discuss the reason for the dif- 
ferent opening angles of the radiation driven outflow 
and the associated different mass loss rates out of the 
pre-stellar molecular cores (Sect. 5.1) and the reasons 
for a longer epoch of sustained disk accretion when gas 
opacities are properly taken into account (Sect. 5.2). 

5.1. Envelope evolution 

Within the low-density optically thin bipolar re- 
gions the choice of gas opacities has little dynamical 
effect. The shape and extent of the low-density bipolar 
region, however, depends strongly on the morphology 
of the accretion disk, where - close to the mid-plane of 
the dust-free zone - the proper treatment of gas opac- 
ities is indeed of importance. In the fiducial run at 
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Fig. 5. — Integrated optical depth of the stellar radia- 
tion at an angle of 15° above the disk's midplane from 
the central star to the outer boundary at t — 45 kyr. The 
optical depth is computed using Planck mean gas and 
dust opacities. The solid and dotted lines denote the 
optical depths of the gas and dust phase, respectively, 
in the fiducial simulation run using the gas opacities by 
Helling et al. (2000), whereby the optical depth is in- 
tegrated radially from the central star to the outer com- 
putational boundary. The dashed and dot-dashed lines 
denote the optical depths of the gas and dust phase, re- 
spectively, in the comparison run with a constant gas 
opacity. 

t — 45 kyr, for example, the dust-free region around 
the forming massive star becomes optically thick to the 
stellar radiation at angles of 9 m +18° above and below 



the disk's midplane. At this point in time, the fiducial 
and the comparison runs differ in their optical depths 
to the stellar radiation (see Fig. 5), resulting in differ- 
ences of the corresponding pre-stellar cores' mass loss 
and their disks' evolution. The gas in the sink cell itself 
is assumed to be optically thin; for further information 
please see the appendix for a comparison to the con- 
vergence simulation run using a smaller central sink 
cell. 

Because the gas is optically thick to stellar radia- 
tion 15° above and below the disk's midplane, dust 
grains persist, where otherwise (without the optically 
thick gas) they would sublimate due to heating by di- 
rect stellar radiation. The radius of dust sublimation 
for an angle of 15° above the disk's midplane at this 
point in evolution decreases from more than 200 AU in 
the comparison run to less than 20 AU in the fiducial 
run. The net effect is to increase the optical depth of 
the innermost 200 AU along the line of sight by more 
than two orders of magnitude. Although it is difficult 
to define precisely the "opening angle" of the outflow 
cavity - it varies both spatially and temporally - we 
note that the large differences of optical depth in the 
inner parts of the envelope lead to an apparent smaller 
opening angle of the outflow in the fiducial simulation 
run (see Fig. 6), which results in a smaller mass loss 
rate of the pre-stellar core. 

The opening angles of the outflow in the fiducial 
run with proper gas opacities, are more consistent with 
observations, whereby molecular outflows tend to have 
opening angles of 90° within 50 AU of the star (Arce 
et al. 2007). 

5.2. Disk evolution 

Disk accretion occurs over a longer period of time 
in the fiducial run in contrast to the comparison run 
with lower (constant) gas opacities. The reason for this 
is two-fold: The so-called "flashlight effect" is greatly 
enhanced by the added opacity contribution of the in- 
ner dust-free gas with its corresponding shift of dust 
sublimation to smaller radii in the shielded disk re- 
gions. This prevents the disk from being torn apart 
by radiative forces for more than t > 200 kyr, i.e. most 
probably for its entire lifetime. In addition, the accre- 
tion disk itself continues to accrete material from the 
surrounding molecular core for a longer period of time, 
because the shielding reduces the ejection of molecu- 
lar gas by radiative forces. Because of the reduced ac- 
cretion onto the disk in the comparison run, it is less 
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Fig. 6. — Gas density distribution of the disk feeding envelope region at t = 63 kyr in the comparison run (upper 
panel) and the fiducial simulation run (lower panel). The snapshot in time is chosen shortly prior to the dispersion of 
the accretion disk in the comparison run. 
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Fig. 7. — Radial gas density profiles at the disks' midplanes at different evolution times as indicated. The solid lines 
represent the fiducial run; the dashed lines denote the comparison run using a constant gas opacity. 
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Fig. 8. — Radial gas temperature profiles at the disks' midplanes at different evolution times as indicated. The solid 
lines represent the fiducial run; the dashed lines denote the comparison run using a constant gas opacity. 
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massive than the fiducial accretion disk. As the den- 
sity in the disk's midplane drops, the ratio of radia- 
tive to gravitational forces increases. Radiative forces 
eventually dominate and drag the remnant disk from 
the central star, and accretion stops (see Figs. 7 and 8 
for a visualization of the density and temperature evo- 
lution in the disks' midplane for both simulation runs, 
respectively). 

5.3. Accretion variability 

The results in this work do not show episodic accre- 
tion analogous to that obtained by Bell & Lin (1994), 
who studied the FU Ori phenomenon using a tempera- 
ture dependent opacity in the disk. However, although 
the opacities used in this study also depend sensitively 
on the temperature, the variations are not as strong as 
in the study by Bell & Lin (1994), who examined the 
thermal instability occurring at the ionization front. 
Moreover, what little variability that we obtained in 
our fiducial run (see Fig. 4) was reduced when we in- 
creased the resolution at the inner disk edge (see dis- 
cussion of the "convergence run" in the appendix). We 
therefore conclude that our accretion variability is nu- 
merical rather than of physical origin. 

6. Observational implications 
6.1. The upper mass limit of stars 

Kuiper et al. (2010a) address the question how mas- 
sive stars can form by disk accretion up to the observed 
stellar mass range in spite of strong radiation pressure 
feedback. In these previous simulations the optically 
thick dusty disk leads to an anisotropy in the thermal 
radiation field (the "flashlight effect") enabling steady 
disk accretion along the radial direction while the ma- 
jority of thermal radiation escapes in the direction per- 
pendicular to the disk's midplane. In this study we 
additionally take into account the dust-free optically 
thick gaseous disks around massive stars interior to the 
dusty disks. This leads to an amplification of the flash- 
light effect. For the collapse initial conditions studied 
here (M core = 100 Mq), the final mass of the star in the 
fiducial run (taking into account the gas opacities) is 
20% higher than in the comparison run (using a con- 
stant gas opacity). Hence, we conclude that the radia- 
tion pressure feedback can efficiently be circumvented 
by the flashlight effect and radiation pressure feedback 
alone does not explain the upper mass limit of stars 
M, < 150 M . 



6.2. Dust sublimation front 

Kraus et al. (2010) report an observation of dust as 
close as 10 AU around a massive star of about 20 M 
(Grave & Kumar 2009). Such small dust sublimation 
radii require efficient shielding of the high-mass stellar 
irradiation. As depicted in our simulation analysis in 
Fig. 5, this shielding can be provided by an optically 
thick gaseous disk. The location of the dust sublima- 
tion front will furthermore be affected by e.g. chemical 
evolution, dust-gas decoupling, magnetic fields, and 
ionization, which we have not included in this initial 
study. 

6.3. Disk lifetimes 

Taking into account the absorption by gas in the 
fiducial run, the feeding of the disk by the large scale 
envelope allows for disk accretion epochs much longer 
than computed in the comparison run using a constant 
gas opacity (at least by a factor of three to four, see 
Fig. 4). We infer from these results that radiation pres- 
sure feedback alone seems to be ineffective in limiting 
the lifetime of accretion disks around massive stars to 
the observational estimates of 1 - 2 x 10 s yr (see e.g. 
Cesaroni et al. 2007). 

On the other hand, the lifetime of accretion disks 
around massive stars might additionally be affected 
by ionization, photo-evaporation, and fragmentation. 
Taking into account the optically thick gas around 
massive stars yields cooler accretion disks, see Fig. 8. 
A cooler disk is more prone to gravitational instabili- 
ties and more likely to fragment. However, fragmenta- 
tion is not studied in the axially symmetric simulations 
performed here. 

6.4. Outflow morphologies 

Although the majority of outflows observed are 
likely to be dominantly driven by magnetic forces, ra- 
diative acceleration may play a role for higher lumi- 
nosity stars in shaping the outflow morphology. In this 
study it is shown that the opening angle of the outflow 
depends on the spatial scales at which it is observed 
and its point of evolution in time. In general, the open- 
ing angles in the fiducial run are smaller than in the 
comparison run and are more consistent with observa- 
tions, whereby molecular outflows tend to have open- 
ing angles of 90° within 50 AU of the star (Arce et al. 
2007). 
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7. Limitations & Outlook 

In the low-density bipolar regions the choice of gas 
opacity is irrelevant because they are extremely opti- 
cally thin. In the dust-free parts of the circumstellar 
disk, however, the correct treatment of gas opacities is 
important; it influences both the large scale shielding 
of the outer regions of the molecular core and the flash- 
light effect. The former allows accretion onto the cir- 
cumstellar disk and the latter allows radial flow within 
the disk. 

This study emphasizes the importance of the inner- 
most disk morphology for the radiative feedback effi- 
ciency. We consider it to be a first step in a broader 
investigation, because many details of the gaseous 
environment and disk morphology close to accreting 
(proto)stars are poorly known. Other effects, such 
as ionization, fragmentation and complex magnetic 
fields, will also influence the local environment. In 
fact, a more comprehensive scan of the broad parame- 
ter space - effecting especially the disk morphology 
such as the initial angular momentum distribution - 
will be necessary to further enlighten the quantitative 
importance of the contribution of the gas phase to the 
properties of high-mass star forming regions. 

Admittedly the enhancement of the flashlight effect 
detected in this study may be even higher, if the inner 
sink cell had been smaller and the spatial resolution in 
the inner regions had been greater, i.e. if more of the 
inner gas disk had been included in the computational 
domain and had been able to contribute to the optical 
depth of the irradiated accretion disk. Based on the 
results of the convergence run with a smaller sink cell 
radius of R m \ n = 5 AU (see appendix), however, we 
are confident that this additional enhancement would 
not significantly alter our conclusions. 

Recent progress of understanding the flashlight ef- 
fect has been achieved by considering more sophis- 
ticated algorithms for radiation transport in the pre- 
stellar core collapse simulations. In particular, the ex- 
tension toward frequency-dependent stellar irradiation 
using a ray-tracing method revealed important effects 
such as the stability of radiation-pressure-dominated 
outflow cavities (Kuiper et al. 2012) and the valida- 
tion of the flashlight effect (Kuiper et al. 2010a). Hav- 
ing this in mind, we could expect to achieve a more 
detailed insight into the radiative feedback of massive 
stars by extending this study to a frequency-dependent 
treatment of gas opacities. Unfortunately, however, be- 
cause of the contributions of molecular line absorption 



and the complex time-dependent molecular chemistry 
involved, this is not straight-forward. 
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A. Convergence run 

To study the dependence of the revealed effects on the size of the sink cell at the inner computational boundary, we 
repeated the fiducial collapse simulation using a sink cell radius of only R m j„ = 5 AU. The gas in the region between 
r = 5 and 10 AU contributes to the total optical depth of the inner gas disk and therefore could alter the results obtained 
in this study. Fig. 9 compares the most relevant physical properties of the collapsing system to the fiducial run. Indeed, 
as expected, the increase of the total optical depth to stellar illumination further enhances the flashlight effect and the 
shielding of the large scale mass reservoir in the molecular core. In addition, the accretion rate into the inner sink cell 
shows less (numerical) variability due to the improved polar resolution at the inner computational boundary. 
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Fig. 9. — Results of the convergence study using a smaller sink cell radius, which allows a larger part of the innermost 
gas disk to be included within the computational domain. Here, the mass loss of the pre-stellar molecular core by 
radiative forces (upper panel), the mass growth of the central massive star (middle panel), and the accretion rate onto 
the central star (lower panel) as a function of time are given. The results of the fiducial run using a sink cell size of 
^min = 10 AU are shown in black. The results of the convergence run using a sink cell size of R m [„ = 5 AU are shown 
in red. 
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